Entropy of polydisperse chains: solution on the Husimi lattice 
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We consider the entropy of polydisperse chains placed on a lattice. In particular, we study a 
model for equilibrium polymerization, where the polydispersivity is determined by two activities, 
for internal and endpoint monomers of a chain. We solve the problem exactly on a Husimi lattice 
built with squares and with arbitrary coordination number, obtaining an expression for the entropy 
as a function of the density of monomers and mean molecular weight of the chains. We compare 
this entropy with the one for the monodisperse case, and find that the excess of entropy due to 
polydispersivity is identical to the one obtained for the one-dimensional case. Finally, we obtain a 
distribution of molecular weights with a rather complex behavior, but which becomes exponential 
for very large mean molecular weight of the chains, as required by scaling properties which should 
apply in this limit. 

PACS numbers: 65.50.+m,05.20.-y 



I. INTRODUCTION 

The statistical mechanics of chains placed on lattices 
has a long history, among the pioneering studies we may 
mention the study of dimers on two-dimensional lattices 
as a model for adsorption of diatomic molecules on a sur- 
face In the simplest version of the model, the dimers 
are chains with two monomers which occupy edges be- 
tween first neighbor sites of the lattice and only excluded 
volume interactions are considered, so that the problem 
is athermal. In the sixties, the entropy of dimers placed 
on two-dimensional lattices was calculated exactly in the 
particular case in which the lattice is fully occupied p|. 
The extension of this calculation to three-dimensional 
lattices or to situations where the lattice is not fully occu- 
pied are still open problems. An interesting recent con- 
tribution in the dimer problem is the extension of the 
solution to the situation in which one site on the border 
of the lattice is not occupied Q, that is, the problem of 
dimers and one monomer. In another generalization of 
the dimer problem a dimer is associated to different ener- 
gies if it is placed on different edges of an anisotropic lat- 
tice [1]. These models lead to unusual phase transitions 
and may be applied to the study of phase transitions in 
some ferroelectric crystals (see [5[ for an example). 

Another generalization is to consider linear chains with 
a number of monomers (molecular weight) M. If the 
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chains are totally flexible, all allowed configurations will 
have the same energy and again the model is ather- 
mal. One fundamental equation, in the thermodynami- 
cal sense, for this gas of M-mers is the entropy per site 
as a function of the fraction of lattice sites occupied by 
monomers, which may be defined as: 

s M (p) = Jim ln t r ^ = P V / M > M > V ^ W 

where T(N P , M; V) is the number of ways to place N p 
linear chains with M monomers each on the lattice with 
V sites. In two dimensions this problem has bee stud- 
ied using series expansions @|, exact solutions on Bethe 
and Husimi lattices Q and transfer matrix solutions of 
the model defined on finite strips followed by finite size 
scaling extrapolations to the two-dimensional limit Q. 

Here we study a version of this model where the molec- 
ular weight of the chains is not fixed. In other words, in- 
stead of the monodisperse set of chains considered above, 
polydisperse chains, with different numbers of monomers, 
are allowed. The particular model we consider is defined 
in the grand-canonical ensemble and was used to model 
equilibrium polymerization. It was applied by Wheeler 
and co-workers to the polymerization transition of sulfur 
The statistical weight of a configuration of chains 
on the lattice is determined by the activities Zi and z e 
for internal and endpoint monomers of chains, respec- 
tively. The particular case z; = corresponds to the 
dimer model (no internal monomers allowed) , and in the 
limit z e — > only very long chains are present. This 
second limit was considered in detail in Q , since a poly- 
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merization transition occurs in this limit. Actually, the 
model may be mapped into the n-vector model of mag- 
netism in the limit n — > and the magnetic field in the 
magnetic model is associated to the activity of endpoint 
monomers z e [9(. 

The solution of statistical mechanical models on hier- 
archical lattices such as the Bethe lattice is a useful 
method to estimate the thermodynamic behavior of these 
models on real lattices [ll| > and the equilibrium polymer- 
ization model was solved on this lattice, both with and 
without dilution fl2l ]. However, this work concentrated 
on the polymer limit of the model, where the polymer- 
ization transition occurs. Recently, we studied the Bethe 
lattice solution of the model in the general case [l3| . Be- 
sides obtaining the entropy as a function of the density 
and the mean molecular weight M , we also studied the 
distribution of molecular weights, which was found to 
be exponential, similar to the distribution for the one- 
dimensional case Here we present a solution of the 
same model on a Husimi tree built with squares. On the 
Bethe lattice, which corresponds to the core of a Cay- 
ley tree, no closed loops are present, and this is one of 
the relevant differences between this lattice and regular 
lattices. On the Husimi tree we considered, small loops 
with four edges are present, and therefore we may ex- 
pect that the solution of models on this lattice should 
be closer to the ones for hypercubic lattices than the one 
on a Bethe lattice with the same coordination number, 
since the elementary plaquette of both lattices will be the 
same. Another motivation for the present calculation is 
that the exponential molecular weight distribution which 
was always found on the Bethe lattice solution may be an 
artificial result for this lattice. Actually, since the critical 
exponents for such hierarchical lattices are classical, an 
exponential distribution is expected in the critical con- 
dition at the polymer limit, since the classical value for 
the exponent 7 is equal to unity fl5j , but deviations from 
this distributions are possible in the general case. 

In section [TT] we define the model more precisely and 
obtain the entropy for a general Husimi lattice built with 
squares, as well as the distribution of molecular weights. 
Final comments and discussions may be found in section 
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II. DEFINITION OF THE MODEL AND 
SOLUTION ON THE HUSIMI LATTICE 

We consider a grand-canonical model of chains placed 
on a lattice. Each chain is composed by two endpoint 
monomers, placed on the lattice sites, and 0,1,2,... in- 
ternal monomers, corresponding to linear self- and mutu- 
ally avoiding walks on the lattice. The statistical weight 
of a particular configuration of chains will be z^"z i i % 
where N e and iVj are the numbers of endpoint and inter- 
nal monomers in the configuration, respectively, while z e 
and Zi are the activities of the two types of monomers. 
The partition function of this model on a lattice with V 



sites may be written as: 

E(z e ,Zi;V) 



(2) 



where the sum over all configurations of the chains on the 
lattice. In Fig. Q]a possible configuration of the chains is 
shown. The density of endpoint monomers is p e = N e /V, 




FIG. 1: A configuration of chains placed on a Husimi tree 
with ramification of squares a = 1. Internal monomers are 
represented as white circles and endpoint monomers are black 
circles. The statistical weight of this configuration with 3 



the density of internal monomers is pi = Ni/V, and the 
total density of monomers is p = p e + pi. The densities 
may be obtained from the partition function as follows: 



and 



z % 

Pi = — -5— , 

a dzi 



z e 

5 dz K 



(3) 



(4) 



Let us now consider the model defined on a Husimi tree 
(or cactus), which is a Cayley tree built with squares. 
The ramification of squares will be equal to er, so that 
the coordination number of the tree is q = 2(a + 1). In 
Fig. [T] a tree with tree generations of squares is shown. 
The tree has a hierarchical structure which allows many 
statistical mechanical models to be solved in quite simple 
ways, similar to the ones used on Cayley trees and Bethe 
lattices To solve the model on a Husimi tree, it is 
convenient to consider rooted subtrees, with a square at 
the root connected to 3cr other subtrees with one genera- 
tion less. The configuration of the two bonds incident on 
the root site of the subtrees is fixed, and we define par- 
tial partition functions (ppf 's) on these subtrees for fixed 
root configurations. The possible root configurations are 
depicted in Fig. [2] so that partial partition function go 
corresponds to a root configuration without any bond 
coming from above, <?;, i — 1,2,3,... are related to root 
configurations with a single bond coming from above, 
connected to i monomers. Finally, the partial partition 
function h corresponds to a root configuration with two 
incoming bonds. The root configurations with one in- 
coming bond have been split into the cases where the 
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bond is already attached to i monomers since we later 
will be interested in obtaining the distribution of sizes 
of the chains. It will be useful to define the sum of all 
partial partition functions with a single incoming bond: 



* = £•■ 



(5) 



g Q 8j.82.g3,... h 

FIG. 2: The root configurations of the subtrees 

Now we obtain expressions for the partial partition 
functions of subtrees with n + 1 generations of squares in 
term of the ones of subtrees with n generations of squares. 
This may be done considering the operation of attaching 
3 groups of a smaller subtrees to a new root square. As 
an example, the recursion relation for a subtree without 
any bond incident on the root site will be: 



fJo 



I a 1 ; 1 °'( Cr 1) cr-2 7 2 

(g +az e g k H Zi g k 



aZihy 



{z e g% + az.kg^ 1 ) 2 [2(g$ + az.g^k + 

-.:;</,; 2 k 2 + azih) + Zig%]. (6) 
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The first term in the sum corresponds to a configura- 
tion of the root square without any polymer bond on its 
edges. Similar recursion relations may be obtained for 
the other ppf's, but since they are rather long they will 
not be written down here. When these recursion relations 
are iterated, we obtain the ppf's of a subtree with addi- 
tional generations of sites. As expected, these partition 
functions diverge in the thermodynamic limit. 

Let us now define the ratios of partial partition func- 
tions: 
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i = 1,2, 



i=l 



(7a) 
(7b) 

(7c) 



As seen in the explicit recursion relation presented above, 
some linear combinations of partial partition functions 
appear repeatedly, so it is convenient to define them here: 



9lA = gZ 



a{<7 — 1) r, 

1 + az e K + — '- Zi K 2 



+crziH] , 
b = g%B = g Q T (z e + az i K), 



(9) 



The recursion relations for the ratios of partial partition 



functions are: 

2z e (A 2 + B 2 ) 



G[ = 



G' = 



G' = 



G . 



H' = 



D 

2z t [a(A 2 +B 2 )d +z e A] 



D 

2z l [a(A 2 + B 2 )G 2 + uZjAGx + ZjZ e ] 
D 

2a Zl [(A 2 + B 2 )G^ + Zl AG t _ 2 + z 2 G^ 3 ] 



(10a) 
(10b) 



(10c) 



D 



(i > 4), 
(A + 2z t )B 2 
D ' 



(lOd) 
(lOe) 



where D = A 3 + B 2 (2A + zt). Summing the recursion 
relations for the Gi leads to: 



K' = J2 G * 



2B(A 2 + B 2 + Zi A + zf) 
D 



(11) 



In the thermodynamic limit, the behavior of the model 
will be described by the fixed point values H' = H,K' — 
K of the pair of recursion relations Eqs. (|10e|) and (jlip . 
Due to the non-linearity of these equations, we were not 
able to find closed expressions of the fixed point values 
H(zi, z e ] a) and K(zi, z e ; a) in general, but their numer- 
ical determination is straightforward. Two particular 
cases where the fixed point equations simplify are when 
only dimers are allowed {zi = 0) and the polymer limit 
[z e — > 0). In these limits the recursion relations above 
reduce to results in the literature 0- Once the fixed 
point values of K and H are found, it is easy to deter- 
mine the remaining ratios G\, G2, . . . using the equations 
above. The partition function of the model on the Husimi 
tree may then be obtained considering the operation of 
attaching a + 1 subtrees to the central site: 



9o +1 



\ + {a + l)[z e K + z l H+-z l K 2 



(12) 



We are not interested in the solution of the model on the 
whole Husimi tree, which is expected to show a behavior 
which is quite different from the one on regular lattices, 
and thus concentrate our attention on the behavior at 
the core of the tree (Husimi lattice). Thus, we obtain the 
densities of external and internal monomers at the central 
site, which are easily found considering the contributions 
to the partition function in Eq. (fT2"|) above. The results 
are: 



-K 



Pe = 



^ + z e K + Zi H + % Zl K 2 '■ 



and 



ZiH + § Zi K 2 



jrrj + z e K + ZiH + %ZiK 2 



(13) 



(14) 
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The entropy per site s(p e , pi) is related to the activities 
through the equations of state: 



and 



In ; 



In Zi 



ds 

dp e 



ds_ 

dp t 



(15) 



(16) 



The entropy may then be obtained inverting Eqs. (|13[) 
and (|14[) to obtain z e {p e , pi) and Zi(p e , pi) and then per- 
forming the integration: 

rpe rf>i 
s(p e ,pi) = - lnz e (p,0)dp- In Zi(p e ,p)dp. (17) 
Jo Jo 

The first integral may be performed analytically, since 
in the dimer limit the fixed point equations are exactly 
solvable. The resulting expression is Q: 

s( Pe ,0) = -(l-p e )ln(l-p e )-^ln^ + 

4(I-o-H-T ' 



ill 

2 V2 

— ~ ln(l + 2W - W 2 ), 



p e ) ln(l + W) 



where 



W 



q 

2 Pt 



q 

2p e 



1/2 



(18) 



(19) 



and q = 2(er+l) is the coordination number of the lattice. 
Now the entropy may be found performing the second 
integration in Eq. (jTTJ) numerically. 

We will, however, proceed along a different path, which 
leads to the same results but with a simpler numerical 
calculation. We start with an argument proposed by Gu- 
jrati [ill ] to obtain the bulk grand-canonical free energy 
per site 6^, adapting it to a Husimi lattice built with 
squares [161 ]. This result may be obtained assuming the 
free energy per elementary square to be a function of 
the generation of the plaquette in the tree and that in 
the bulk it will converge to a fixed value in the thermo- 
dynamic limit (l7j ]. A simple derivation for the present 
case, where we suppose that the center of the Husimi tree 
is a site, starts by associating a generation number m to 
each site of the tree, such that m — for the central 
site, m — 1 for the 3(cr + 1) first neighbors of the central 
site and so on, until the surface of the tree is reached for 
rn = M. The number of sites at the surface of the tree 
will be equal to: 



The number of remaining sites in the tree, which are in 
the bulk, will be: 



M-2 ,„ s 

N b (M) = l+3(cr+l) J" (3aY = 1+3(^+1) 1 a) 



M-\ 



- 1 



3=0 



3(7 - 1 



(21) 

Supposing that the free energy of the whole tree may be 
written as a sum of the contributions of its surface and 
the bulk, we have: $(M) = -k B TlnE(M) = N b {M)(j) b + 
N s (M)(j> s , where 4>b and 4> s are the free bulk and surface 
free energies per site, respectively. In the thermodynamic 
limit M — > oo, we may then find that: 

-lnH(M + l) + 3o-lnH(M) = [N b (M + l)-3aN b (M)]tp b 

(22) 

where ip b — ^/(fcsT). The substitution of the expres- 
sions pO|) and ([2Tj) for the numbers of bulk and surface 
sites we will find that: 



1 f Z(M + l) 



(23) 



This expression is a generalization of the one obtained 
using a similar argument by Semerianov and Gujrati for 
(j = l (expression A7 in [16(]). If we now substitute the 
partition function (|12p in the expression above and ob- 
serve that go(M + l)/(go) 3<7 is equal to the fixed point 
value of the denominator of the recursion relations (jlOej) 
D, so that 



D a+1 



[1 + ((7 + \){z e K + Zi H + az^/2)}^- 1 



<Pb = 
,(24) 



where the thermodynamic limit is implicit since the fixed 
point values of D, K, and H are used. 

The entropy per site in the bulk of the tree is given by 
the state equation: 



cm 
dT 



which leads to the expression: 



t— = -<fb - p e lnz e 
Kb 



Pi In Zi 



(25) 



(26) 



Solving equations (1131) and ([14]) for the activities z e and 
Zi, we obtain the equations: 



and 



Pe = z e K (a + 1)(1- Pe- Pi), (27) 



Pl = Zi (a + 1)(H + aK 2 /2){\ - Pe - Pi ). (28) 



N S {M) =3(cr + l)(3cr) 



M-l 



(20) 



To obtain the entropy as a function of the densities of 
internal monomers and endpoint monomers, we numer- 
ically solve the set of four equations (|2"7|) . (|2"5|). (|10cl) . 
and (jlip (making H' = H and K' = K in the last two 
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equations), so that we obtain z e , Zi, H, and K as func- 
tions of the densities, and then the entropy per site in 
the bulk using expression (|26[) . In general, this numer- 
ical procedure converges to the fixed point and leads to 
accurate results, but there are problems when the full lat- 
tice limit is approached, as expected, since the activities 
diverge in this limit. Therefore, to compute the entropy 
in the full lattice limit we use a separate procedure we 
will describe below. We thus take the limit of diverg- 
ing fugacities keeping z — z e /zi — exp[(/i e — p^ / (ksT)] 
fixed. In this limit, the variables defined in expressions 
© and © diverge and it is convenient to define the new 
variables 



D 



A = azK+ a[a l) K 2 + aH, (29a) 
Zi 2 

B 



z + aK. 



^=A 3 0O +B 3 O0 (2A 0O + l), 



and the recursion relations may be rewritten as 



H' = 



K' 



BUAoo + 2) 



Do 

Ml 



'OO 

2B n 



A n 



1) 



D c 



(29b) 
(29c) 

(30) 
(31) 



The density of endpoint monomers will be given by: 

zK 

p e = 



zK + H + <jK 2 /2' 



(32) 



and we notice that in the limit z — > oo, which corresponds 
to the lattice fully covered with dimers, we find that p e — > 
1, while for z — > we have p e — > 0, as expected. The 
density of internal monomers will be pi — 1 — p e , and 
the mean number of monomers per chain is M = 2/ p e . 
Substitution of equation f|32[) in this last expression leads 
to: 



H = ^[z{M-2)-aK]. 



(33) 



This last expression allows us to eliminate H from the 
fixed point equations associated to the recursion relations 
(I10e|) and ([TT1) . so that the fixed point equations reduce 
to: 



A", 



2 

oo ' 



[z(M - 2) - aK]Doo = (A^ + 2)B, 
KDoc - 2B 00 [A 00 (1 + Aoc) +^ + 1] 



(34) 
(35) 



For fixed values of the mean molecular weight of the 
chains M, these equations may numerically be solved 
for K and z, using also the equations (|29|) . Once the 
fixed point is found, we may obtain the entropy noticing 
that expression (P2B1) in the full lattice limit reduces to 
s/kB = —ifb — Pe In z — In Zi, and substitution of the bulk 
free energy in this limit allows us to write the entropy as: 



— In — — 

4 [(*■ 



D <T+1 



l)MzK/2} 3 < 7 - 



M 



■ In z. 



(36) 




FIG. 3: Entropy as a function of the density of monomers 
p — p e + Pi for polydisperse chains on the Husimi lattice. In 
the region of small densities, the curves shown correspond to 
mean molecular weights M = 3, 6, and 10 in downward order. 
The results are for a = 1 (q = 4). 




FIG. 4: Entropy as a function of the mean molecular weight 
at full coverage for a Husimi lattice with a = 1 (q = 4). 



On the Bethe lattice, it is possible to obtain an the 
entropy explicitly using the expression for the bulk free 
energy, and we show this calculation in some detail in the 
Appendix. 

Results for the entropy as a function of the total den- 
sity of sites, for fixed values of the mean molecular weight 
M = 2{p e + Pi)/p e are displayed in figure [3] As was 
already noticed in the results on the Bethe lattice, the 
entropy is a increasing function for small values of the 
density and goes through a maximum at intermediate 
densities. We notice that at low densities the entropy is 
a decreasing function of the mean molecular weight M, 
but this does not hold for larger densities. 

It is interesting to compare the present results with 
the ones for the Bethe lattice [l3j]. In general, as was 
also found in the monodisperse case Q, both results are 
quite close. In figure [5] the relative differences between 
both entropies are given as functions of the density for 
some values of M, in general these differences are in a 
range below 1%. At low densities, the entropy on the 
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p 



FIG. 5: Relative differences between the entropies of poly- 
disperse chains on Bethe (sb) and Husimi (sh) lattices as 
functions of the density for lattices with q — 4. From the up- 
per to the lower curve, the mean molecular weights are equal 
to 10, 6, and 3. 




P 



FIG. 6: Relative differences between the entropies of poly- 
disperse chains on Bethe (ss) and Husimi (sh) lattices as 
functions of the density for M = 6. From the upper to the 
lower curve, the coordination numbers are equal to 4, 6, 8, 
10, and 12. Numerical round-off errors are already visible at 
higher coordination numbers, particularly at low densities. 



Bethe lattice is larger than the one on the Husimi lattice, 
but the opposite may happen at higher densities. The 
difference between the results on both lattices reduces 
as the coordination number q is increased, as expected, 
since in the limit q —y oo both entropies are equal to the 
simple mean field estimate, this is apparent in the plots 
displayed in figure HJ 

In the table U we present data for the entropy at full 
coverage (p = 1). As mentioned before for the general 
case, the differences, both in the monodisperse and the 
polydisperse systems, between the results of the Bethe 
lattice and of the Husimi lattice solutions are quite small, 
and there are indications that these values are still rather 
far away from the results for regular lattices with the 
same coordination number. For dimers (M = 2), where 
there is no polydispersivity in our model, the exact value 



M Bethe (m) Husimi (m) Bethe (p) Husimi (p) 



2 


0, 26162 


0, 26743 


0, 26162 


0, 26743 


3 


0, 42284 


0,41295 


0, 88493 


0, 88650 


4 


0,48166 


0,48951 


0, 95904 


0, 95969 


5 


0, 50669 


0, 50888 


0, 95656 


0, 95673 


6 


0, 51349 


0,51265 


0, 93473 


0, 93467 


7 


0,52217 


0, 52284 


0, 90836 


0, 90820 


oo 


0, 4055 


0, 4090 


0, 4055 


0, 4090 



TABLE I: Entropy at full coverage for lattices with q = 4. 
Results for monodisperse (m) and polydisperse (p) chains are 
shown. Data for the monodisperse case are from reference 
[3 and the case of polydisperse chains on the Bethe lattice is 
discussed in 

on the square lattice is known Q (s2 = G/tt 0.29156, 
G is Catalan's constant). In this case, the value on the 
Bethe lattice is about 10% below the exact value, while on 
the Husimi lattice the difference reduces to roughly 8%. 
The entropy for monodisperse chains with M > 2 calcu- 
lated on Bethe and Husimi lattices is always larger than 
the estimates obtained from transfer matrix calculations 
for the square lattice [§[, and the relative differences are 
smaller than the ones for dimers. In the polymer limit 
M — y co, the difference between the entropies for the 
poly- and monodisperse cases vanishes again. 

We notice that the entropy at full coverage is not a 
monotonic function of the mean molecular weight, show- 
ing a maximum around M « 4.37, as may be seen in fig- 
ure|U As expected, on the Husimi lattice we find that the 
entropy for polydisperse chains is always larger or equal 
to the one for monodisperse chains with M monomers for 
the same density of monomers p and finite M = M > 2, 
obtained in [7J, where equality holds only for vanishing 
density. 

On the Bethe lattice, it was possible to obtain the 
entropy for polydisperse chains analytically, and it was 
found that the contribution of polydispersivity to the en- 
tropy, 

As M (p) = s M (p) - s M (p), (37) 

with M — M, is linear in p and independent of q f]~3j j : 
As m ,b = p[(M - 1) ln(M - 1) - (M - 2) ln(M - 2)]/M. 
The contribution to the entropy of polydispersivity in the 
solution of the model on the Husimi lattice (Asm,b does 
not show such a simple behavior. In figure [7] we show 
the difference between the results on the Bethe and on 
the Husimi lattices, for different values of the molecular 
weight M. We notice that in general, as the molecu- 
lar weight M is increased, the Husimi lattice results are 
closer to the ones found on the Bethe lattice. This may 
be understood if we remember that since no loops are 
present on the Bethe lattice and loops of four edges only 
may be found on the Husimi lattice we considered, a rea- 
sonable effect of these closed paths should be expected 
for chains of length close to four, but these effects should 
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FIG. 7: Difference between the contribution of polydispersiv- 
ity on the Bethe and the Husimi lattice solutions as a function 
of the density p. The full curve is for M = 3, the dashed one 
for M = 4 and the dotted one for M = 5. All results are for 
a = 1 (q = 4). 




FIG. 8: Difference between the contribution of polydispersiv- 
ity on the Bethe and the Husimi lattice solutions as a function 
of the density p. The full curve is for a = 1, the dashed one for 
a = 2 and the dotted one for a = 3 (q = 6, 4, 8, respectively). 
All curves are for M = 4. 



become smaller as the chains grow. In opposition to what 
is found on the Bethe lattice, the contribution of poly- 
dispersivity to the entropy on the Husimi lattice solution 
changes as the ramifications of squares a varies. This 
may be seen in figure [SJ In general, we find that the 
Husimi lattice results become closer to the ones found 
on the Bethe lattice as a, and therefore the coordination 
number q, grow. This is expected, since both solutions 
should become equal to the simple mean field solution in 
the limit q — > oo. However, as is apparent in the curves 
in figure [SJ the convergence is not monotonic. 

To calculate the probability tm to find a chain with M 
monomers among all chains, we notice that: 



G 



M-l 



K 



(38) 



and this probability can be numerically evaluated, since 
we may obtain the functions G and K, whose recursion 



relations are given in Eqs. (jlOel) and (ITTT) . at the fixed 
point. We recall that on the Bethe lattice these ratios are 
given by B B =r M +i/r M = (M-2)/(M-l), and there- 
fore are independent of p and q [l3[ . As may be noted in 
Fig. [9l where the ratios Bh(M) = Tm+x/^M are plotted 
as functions of M, for small values of M compared to 4 
the distribution of molecular weights is not exponential, 
but as M 3> 4 an exponential behavior is apparent. The 
value 4 corresponds to the size of the only possible closed 
path on the lattice we considered. Again, we notice that 
the deviations from the exponential behavior are larger 
for smaller values of q, as q increases the results approach 
an exponential behavior for all values of M , as found on 
the Bethe lattice. The asymptotic value of the ratio for 
large values of M may be found by assuming an exponen- 
tial decay of the probabilities for the fixed point values 
of the ratios of partial partition functions Gi for i > 3 in 
Eqs. (|10e|) . which leads to the following equation for the 
limiting ratio B H = lim M ^oo r M +i/r M - 



B) B 2 H + Zl A B 



H 



= 0, (39) 



and the horizontal line in Fig. [S] was obtained solving 
this equation. 




FIG. 9: Ratios of probabilities to find chains with successive 
molecular weights M on the Husimi lattice Bh(M) as func- 
tions of the molecular weight M for M — 3. Data shown 
are for a = 1 and p = 0.3019045. The dashed line cor- 
responds to the result on the Bethe lattice, which for this 
case is Bb = 1/2. The full line is the asymptotic value 
B H =lim A /^ 00 B H (M). 

The asymptotic ratio of the solution on the Husimi 
lattice, Bh-, in general will be a function of M, a and p. 
In figure [TU] this is apparent, and again we notice that as 
the coordination number of the lattice become larger, the 
result on the Husimi lattice approaches the one found on 
the Bethe lattice. This may be seen analytically in the 
limit of vanishing activities. Expanding the recursion 
relations up to the lowest nonzero order of z e and z,, we 
get: 



H' 
K 1 



z\ + 2(7Z e Z l K + ((TZrK) 2 , 



2z e + laZiK. 



(40a) 
(40b) 
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It is then easy to find the fixed point values up to lowest 
non- vanishing order K* = 2z e and H* = z 2 el and the 
corresponding values for the mean molecular weight M = 
2 + (1 + 2a)zi and density of monomers p = 2(a + l)z 2 . 
The approximate solution of the equation (f39| for the 
decay exponent in the limit of small activities is /3h — 
2a Zi, which may be combined with the expression for M 
above leading to: 

P B = j^(M-2) t (41) 

where the dependence of the exponent with a in the limit 
considered is clear. In general, the exponent will also 
be a function of p, but, as may be seen above, in the 
limit we consider here this contribution is of higher order. 
The same expansion may be done for the Bethe lattice 
solution [Til, and leads to the result /3b = a'zi, where a' 
is the ramification of this lattice. If we compare Bethe 
and Husimi lattices with the same coordination numbers, 
we have a' = 2a + 1, and as expected (3h -> as 
a — > oo. 
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FIG. 10: Asymptotic coefficient Bh as a function of the den- 
sity p, calculated for M = 6 



III. FINAL COMMENTS AND DISCUSSIONS 

We have studied an athermal model of flexible chains 
with excluded volume interactions only on a Husimi lat- 
tice built with squares, with an arbitrary value for the 
ramification of squares a, so that the coordination num- 
ber of the lattice is q = 2(a + 1). The chains are lin- 
ear and composed by a set of monomers, so that con- 
secutive monomers occupy first neighbor sites of the lat- 
tice. They are polydisperse, and the distribution of sizes 
is determined in an annealed way by two parameters: 
the activity z e of an endpoint monomer (linked to one 
other monomer of the chain only) and z$ of an internal 
monomer, which is linked to two other monomer of its 
chain. 

The entropy as a function of the fraction of sites occu- 
pied by monomers and the mean number of monomers in 



each chain (mean molecular weight of the chains) is a fun- 
damental equation of the system in the thermodynamic 
sense, and we have obtained this function in general on 
the Husimi lattice, although, at variance to what was 
done in the solution of this problem on the Bethe lattice 
fl3j , we were not able to derive a closed form for it. We 
found that the differences between results on Bethe and 
Husimi lattices with the same coordination numbers s 
rather small, below 1%, as was also the case for monodis- 
perse chains Q • It is also possible that the results on the 
Husimi lattice for q — 4 are still not very close to the ones 
on the square lattice, although, to our knowledge, no re- 
sults for the entropy of this model on regular lattices are 
available in the literature. For the monodisperse rather 
precise estimates were obtained using the transfer matrix 
approach [8| , and we are presently extending these results 
for the polydisperse model on the square lattice. It is ex- 
pected that the results on the Husimi lattice will be closer 
to the ones on regular lattices at higher dimensions, but 
again there are no estimates of the entropy of the model 
available in the literature for three-dimensional lattices, 
for example. 

Another point which should be further investigated is 
the distribution of molecular weights of the chains. The 
simple exponential distribution which was found on the 
Bethe lattice [l3[ is no longer valid on the Husimi lattice, 
so we expect that on regular lattices a different size dis- 
tribution will be found too. It may be that this question 
could be addressed on the square lattice with transfer 
matrix techniques, we are now investigating this possi- 
bility. 
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Appendix A: Free energy of the model on Bethe 
lattice 

On a Bethe lattice with arbitrary coordination number 
q, the entropy of chains whose polydispersivity is deter- 
mined by different activities for internal and endpoint 
monomers may be calculated exactly. Here we show that 
this result, which was originally obtained integrating the 
expressions for monomer densities, may also be found 
more directly by derivation of the bulk free energy calcu- 
lated using Gujrati's prescription. Recalling the discus- 
sion of the problem in (l3j . we define a partial partition 
functions go and g\ for subtrees with and without a poly- 
mer bond on the root edge, respectively. The recursion 
relation for the ratio R — g\ /go of these ppf 's is given by 
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expression (13) in this paper: 



Bl 



<JZjR 



<rz e R 



z l R 2 ' 



(Al) 



and recalling the recursion relation for go (expression 
(11)), which is: 



and equation (14) for the partition function of the model 
on the Cayley tree may be written as: 



i—i q i 0—1 

= = 9o + 1 z e% 91 



9(9-1) 



9-2 2 

-ZiSo 9i 



9l 



(A2) 



The free energy of the model on the Bethe lattice cor- 
responds to free energy on the bulk of the Cayley tree, 
and denoting this free energy per site of the tree by (f>b 
we obtain, using Gujrati's ansatz, the result [l8l ]: 



4>b 1 -, S m+ i 

LOh = = in 

V k b T 2 



-9-1 
1 — m 



(A3) 



where m denotes the number of generations in the tree 
and we are interested in the thermodynamic limit m — ¥ 
oo. Using expression (|A2[) for the partition function, we 
get: 



fb = 




1 + qz e R - 



9(9-1) 

2 



9-2 



(A4) 



9o = 9o + vz e g 9i + 



a-l , CrO - 1) (7- 



(A5) 



we reach the following expression for the bulk free energy 
per site: 



1 + az e i? + ^i^i? 2 



(l + g z e i?+ aitil Zi i?2) 



g-2 



(A6) 



Now the entropy is given as a state equation associated to 
the free energy, whose expression is equation (f2l)l) above, 
since the activities may be written as functions of the 
densities in the core of the tree, expression (25) in [l3[ 
for the entropy as a function of these densities is found. 
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